Orthogonal dimer-plaquette chain with alternating Ising-Heisenberg coupling 
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We explore the dimer-plaquette Ising-Heisenberg chain, assembled between plaquette edge also 
known as orthogonal dimer plaquette. Using the local gauge symmetry of this model, we are able 
to map onto a simple spin-1 like Ising and spin-1/2 Heisenberg dimer model with single effective 
ion anisotropy, thereafter this model can be solved using the decoration transformation and transfer 
matrix approach. First we discuss the phase diagram at zero temperature of this model, where 
we find five ground states, one ferromagnetic, one antiferromagnetic, one triplet-triplet disordered 
and one triplet-singlet disordered phase, beside a dimer ferromagnetic-antiferromagnetic phase. It 
is interesting to note that, in the limit of pure Ising model there is a frustrated region. In addition, 
we discuss the thermodynamic properties such as entropy, where we display the residual entropy. 
The specific heat also is discussed, where we observe anomalous peaks close to the first order phase 
transition in the low temperature limit and for weak Heisenberg anisotropic coupling, but for strong 
Heisenberg anisotropic coupling, this anomalous peak disappear. Furthermore, using the nearest 
site correlation function it is possible also to analyze the pairwise thermal entanglement for both 
orthogonal dimers, as well as the we discuss the threshold temperature of the entangled region as a 
function of Hamiltonian parameters. 



I. INTRODUCTION 

Quasi two-dimensional magnets have been attracted 
since 90 decades, such as the quasi-two-dimensional mag- 
net CaV409[l], that has a layered structure where the 
magnetic V*"*" ions have spin 1/2 and form a 1/5- 
depleted square lattice. As well as the polycrystalline 
SrCu2 (603)2 having a two dimensional (2D) orthogo- 
nal network of Cu dimers, this cuprate, provides a 2D 
spin-gap system in which the ground state can be solved 
exactly [2, 3]. These quasi- two-dimensional systems are 
topologically equivalent to the theoretical model pro- 
posed by Shastry and Sutherland [J|. 

Motivated by the above real materials Ivanov and 
Richter[5] proposed the class of one-dimensional Heisen- 
berg spin models (plaquette chains) related to the 
real materialspJ-Q, through numerical and analytical 
results was analyzed the zero temperature magnetic 
properties [5i,i^). While in reference [7l| was discussed the 
sequence of first-order quantum phase transitions in a 
frustrated spin half dimer-plaquette chain. A detailed 
investigation about the first-order quantum phase trair- 
sition of the orthogonal-dimer spin chain also was consid- 
ered by Koga et al.[8!|, as well as the frustration-induced 
phase transitions in the spin-S" orthogonal-dimer chain[9]. 

A more recent investigation was developed by Ohanyan 
and Honecker[10], where they have been discussed 
the magnetothermal properties of the Ising-Heisenberg 
orthogonal-dimer chain with triangular XXZ clusters. 
Furthermore, in the last decade several quasi-one- 
dimensional Ising-Heisenberg model such as diamond 
chain were intensively investigated, mainly the ther- 
modynamic properties aird geometric frustration[ll|- 
[l4| . magneto-caloric effect |15|. as well as thermal 
entanglement |l6l Il7| . among other physical quantities. 
Some other variant of the Isiirg-Heisenberg model also 
were considered, such as Ising-Hubbard model[la| and 
Hubbard model in the quasi- atomic limit fl9| besides 



spinless electrons[20] in diamond chain. 

On the other hand, quantum entanglement is one of the 
most attractive types of correlations that can be shared 
only among quantum systems[21]. lir recent years, many 
efforts have been devoted to characterizing qualitatively 
and quantitatively the entanglement properties of con- 
densed matter systems, which are the natural candidate 
for application in quantum communication and quan- 
tum information. In this sense, it is quite relevant to 
study the entanglement of solid state systems such as spin 
chains [221 . The Heisenberg chain is one of the simplest 
quantum systems, which exhibits the entanglement, due 
to the Heisenberg interaction is not localized in the spin 
system. Several studies have been doire oir the thresh- 
old temperature for the pairwise thermal entanglement 
in the Heisenberg model with a finite number of qubits. 
Thermal entanglement of the isotropic Heisenberg chain 
has been studied in the absence [2J| and in the presence 
of external magnetic field [2J| . 

The outline of this work is as follows. In sec. 2 we 
present the dimer-plaquette Isiirg-Heisenberg chain; it is 
also discussed the zero temperature phase diagram. In 
sec. 3 we present the exact thermodynamic solution of 
the model, thus we discuss the eirtropy aird specific heat. 
In sec. 4 we discuss the thermal entanglement and its 
threshold temperature. Finally in sec. 5, we summarize 
our results and draw our conclusions. 



II. ORTHOGONAL DIMER-PLAQUETTE 
ISING-HEISENBERG CHAIN 



The theoretical investigation of the orthogonal dimer- 
plaquette models is motivated not only from the the- 
oretical point of view, but also from the experimental 
viewpoint. It is worth to remark that the Heisenberg 
orthogonal dimer-plaquette model cannot be solved ex- 
actly at finite temperature. Motivated by the comments 
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Figure 1: (Color online) Schematic representation of dimer- 
plaquette Ising-Heisenberg chain. (Above) Thick line corre- 
spond to Heisenberg coupling, while the thin line corresponds 
to Ising coupling. (Below) Dimer plaquette mapping through 
local gauge symmetry. 



given in the introduction, we consider the orthogonal- 
dinier plaquette chain with Ising-Heisenberg coupling as 
schematically described in figure[T] Therefore, the Hamil- 
tonian for an orthogonal dimer-plaquette chain could be 
expressed by 

N 
i=l 

with 
J'icTa,^, <t,,)a' = J'«.<. + 'T^,,<,) + A'a^^,<, (2) 

where a" ^ are the spin operators also known as Pauli 
matrices (with a = {a;,j/, z}) at plaquette i for particles 
7 = {a, 6,c, d}, for detail see figure [TJ The thick line 
correspond to Heisenberg coupling, while the thin line 
corresponds to Ising coupling. The Ising coupling pa- 
rameter is denoted by Jq, whereas J [J') represents the 
X and y components of Heisenberg coupling and with A 
(A') we mean the anisotropic (z-component) coupling in 
the Heisenberg term for dimer ah and cd respectively. 

To transform this model onto the well known mixed 
spin "Ising "-Heisenberg model, we use the following def- 
inition Sf = o'^i + o-^i- By the use of S" definition we 
obtain the identity: (s°'f == 2 + 2cr^,cr^-. Thereafter, 
we can easily establish the following transformation 






{S[f-{2J' + A'), (3) 



and by S^ we just denote S.^ — Si. Si, so, this matrix is 



given by 



SJ=A 



2 

110 

110 

2 



(4) 



rewritten the Hamiltonian in terms of operators S^ 
and S.i we have an Ising-Heisenberg chain model whose 



Hamiltonian is given hy H — Hq -\- H', with Hq = 
(2J' -I- A')N. Therefore, the transformed Hamiltonian 
reduce to 



H' = - 



(SD' + isUiY 



+ J{(Tc,i, eTd,i)A + Jo {SiCTc.i + 0-d,zSt+l)} 



(5) 



We can observe from eq. @ the matrix has a 2 x 2 block 
matrix; hence, we can diagonalize this block matrix. It is 
interesting to note that the matrix S^ is still diagonal in 
the new basis, due to corresponding 2x2 block matrix be 
null. Then, this means we can simultaneously diagonalize 
both matrices Sf and S^. Recall those Sf and S^ are 
commutative operators. So, the diagonal matrices are 
expressed by 



S = 



8 

8 



8 



and Sf = 



2 


0-2 



(6) 



and whose corresponding eigenvector states are expressed 
as 



'+i>=ii>, 


(7) 


ko)-^(|t> + l;)), 


(8) 


l«o)=^(l+)-l;)), 


(9) 


--i)=|:>- 


(10) 



The effective Ising "spin" in Hamiltonian ([SJ can be 
understood as a composition of one triplet state and one 
singlet state. 

This transformation is possible because, the local 
gauge symmetry is only satisfied by Hamiltonian ([1]) 
when the Ising couplings ac (ad) and bc(bd) are identical. 



A. The zero temperature phase diagram 

In order to analyze the phase diagram of the 
orthogonal-dimer plaquette, we assume the particular 
case J' — J and A' = A, following the parameter used 
in the literature^-Q. 

The dimer-plaquette Ising-Heisenberg chain described 
by the Hamiltonian ([T]) exhibits five ground states. These 
states are expressed as 



N 



\FM)=l[\X), 

i=l 

N 

\AFM)^Y[\+), 



+)^ 



\'~)^ 



(11) 

(12) 



»=i 



One ferromagnetic (FM) phase and another antiferro- 
magnetic (AFM) phase, with corresponding eigenvalues 



given by 

EpM = - 2 J' - 2A' - A - 4Jo, 

Eafm = - 2 J' - 2A' - a + 4Jo. 

The other phase are composed by dimers in triplet- 
triplet (TT) phase and triplet-singlet (TS) phase, 



^ 1 

\TT)^\{\t,),® — ^ 

i—1 ^ 



N 



\TS)^\[ 



ro) 



V2 



(13) 
(14) 



where the corresponding eigenvalues become 

Ett =A - 2 J - 4J', 
Ets =A + 2 J. 

It is worthy to mention that the state given by eas. (jl3p 
and P^ becomes frustrated state only when J = (Ising 
limit). 

Finally, the dimer-antiferromagnetic (DFA) state can 
be expressed by 

N/2 
\DFA) =Y[\+),(g, \r,2,^2)^ (E>\Z)^(E> \7l2,-2), (15) 



with 



and 



i=l 



l'72,-2) 



VTTW 



^ = 



y/4J-^ + J2 + 2 Jo 



J 



(16) 



(17) 



Hereafter the corresponding dimer-antiferromagnetic 
eigenvalue becomes 



Edfa - - 2J' - 2A' + a - 2^4 J2 + J^. 

In order to display the phase diagram at zero tem- 
perature, we plot J versus A, assuming fixed value for 
Jo = 1(-1) FM(AFM) respectively. As shown in fig- 
ure [JJa). For A < — 2, there is a boundary between TS 
and TT state at J = 0, whereas for —2 ^ A ^ 0.5 the 
state TS (TT) is limited by DFA region whose boundary 
is given by A = 2| J| - V^Tl^., when J > (J < 0) 
respectively. While for A ^ 0.5 the states TS (TT) is 
limited by FM(AFM) region A = | J| - 1. Furthermore, 
the boundary between DFA and FM(AFM) region is rep- 
resented by the curve A = \/4 -I- J^ — 2. 

While in figure [2l^b), we observe the phase diagram 
from another viewpoint, A against Jp, assuming fixed 
J = 1. In this picture the TS and TT phase is displayed 
in same region, while AFM and FM states are illustrated 




Figure 2: (a) Phase diagram at zero temperature as a de- 
pendency of J and A for a fixed value of Jo = 1. (b) Phase 
diagram at zero temperature as a dependency of Jo and A for 
a fixed value of J = 1. 



for Jo > (Jo < 0) respectively. The boundary be- 
tween FM(AFM) and TT(TS) is given by A = 1 - | Jo| 
for I Jo I < 2/3. Whereas the DFA region appears for 
I Jol > 2/3, which is limited by TT(TS) region and this 
boundary curve becomes A = 2 — ^4Jq + 1, and finally 
the boundary between FM(AFM) region is described by 

2|Jo|. 



the curve A = y/4J^+l 



III. THE THERMODYNAMICS OF THE 
MODEL 

The method to be used will be the decoration trans- 
formation proposed in early 50 decade by Syozi [25| and 
Fisher [20]. Afterward this approach was the subject of 
study in reference (27|, for the case of multi-spins. Similar 
generalization was developed by Strecka [2^ for the hy- 
brid system (e.g. Ising- Heisenberg). Another interesting 
variant of this approach also was discussed previously, 
where was proposed a direct transformation instead of 
several step by step transformation [29]. Consequently, 
the decoration transformation approach is widely used to 
solve spins models, besides, the decoration transforma- 
tion approach can also be applied to electron coupling 
system such has been applied for the case of spinless 
fermion on diamond structure [20], and as well as for Hub- 
bard model in the quasi-atomic limit [19i]. 

In order to study the thermodynamics of the dimer- 
plaquette Ising-Heisenberg chain, we will use the deco- 
ration transformation proposed in reference 129] together 
with the usual transfer matrix technique |30|. So, let us 
start considering the partition function as follow 



-N 



- ^-0H„ 




,-PH' 



(18) 



where /3 = l/kT , with k being the Boltzmann constant 
and T is the absolute temperature, and assuming H' is 
given by eq. (O . 

The effective model can be solved by the usual transfer 
matrix approach]30]. Inasmuch as the transfer matrix of 



the Hamiltonian ([5]) is reduced to 



/2 


/4 

wo,qx' 

WofiX'^ 
/2 


/2 


Wl,-1 
/2 

wi,oa; 


Wl,0 


wo,o 


Wl,0 


wi -1 


Wl,0 


Wl,l 



with 



wi^i =a;''z'' [^ (/ + y-^) + z-\x^ + x'^)] 



X 

wifi ^x'z' [z {y^ + y^"^) + z^^{x\ + x^"^)] 



{y + y 

wi,_i =a;' z' \2z + z 
wo,o =2z + z^ [x - 



')] 



(19) 



(20) 
(21) 
(22) 
(23) 



where we have introduced the following notations x 



=PJ 



p/SA 



,/3J' 



el^^' and 



y 



JiJa 



Furthermore, we define also the following exponentials 
XI = ef'^/^^+^o and X2 = e^^/^^+^o . 

In order to diagonalize the transfer matrix T, we per- 
form the determinant of the transfer matrix, which is 
expressed as follow 

A (A2 - fliA + ao) (A - Wis + wi,-i)) = 0, (24) 

here, the coefficients of the quadratic term eq. ([M)) is 
given by 



ai =wi,i + wofi{l + x' ) + wi,-i, 
oo = f 1 + x' j [wo,o {wis + "^i.-i) 



(25) 



2w 



l,Ol J 



Therefore, the eigenvalues can be expressed as 
Ao =j 



Ol 



Ai =\ Ui 



ai — 4ao 



Aai 



(26) 



(27) 



\ /2 /2 

At =X Z 



Aa =0. 



2(y 



{x" 



(28) 
(29) 



In the thermodynamic limit, the free energy per unit 
cell is given only by the largest eigenvalue of transfer 
matrix, and for our case, we can easily verify the first 
eigenvalues Ao become always the largest one. Thus, the 
free energy can be expressed by the relation below 



1 



/ = 2J' + A'--ln(Ao 



(30) 



The first term of the free energy is a trivial constant 
energy, obtained during the local gauge transformation, 
which is irrelevant for thermodynamis quantities. Note 
that the free energy is valid for arbitrary value of J' 
and A', although here we consider only a particular case 
J' ~ J and A' = A, following the parameter used in the 
literaturefo-Q. Therefore, any additional properties can 
be obtained straightforwardly from eg. pop . 
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Figure 3: Density plot of entropy, darkest region corresponds 
to higher entropy, (a) Entropy as a dependency of J and A 
for a fixed value of Jo = 1. (b) Entropy as a dependency of 
Jo and A for a fixed value of J = 1. 



A. Entropy and specific heath 

In order to accomplish with our discussion concerning 
to the thermodynamics properties, let us illustrate the 
entropy {S = —df/dT) in fig. [3)Ja) for the low tempera- 
ture limit as a function of A and J, whereas, in figEl^b) 
the entropy is illustrated as a function of A and Jq. Dark 
region corresponds to higher entropy, while the white re- 
gion corresponds to zero entropy. However, as soon as 
the temperature decreases the entropy leads to zero, un- 
less, for J = and A ^ — 2 (pure Ising chain) the dark 
line will remain when temperature decreases, leading to 
a residual entropy Sq = 21n(2), that illustrates the ge- 
ometric frustration region of orthogonal plaquette Ising 
chain. 

Finally, we discuss also another interesting ther- 
modynamic quantity called the specific heat (C = 
—Td^f/dT^). In figure 21 we plot the specific heat as 
a function of the temperature T for a fixed coupling pa- 
rameter Jo = 1.0 and a fixed value of J = 1.0. By a solid 
line is represented the specific heat curve when A — 0, 0.5 
and 1.0, while, for dashed line we illustrate the specific 
heat for A — —0.5, and —1.0. For anisotropic parameter 
A = ±0.5 and 0, the plot illustrates a small anoma- 
lous peak due to the zero temperature phase transition 
infiuence, although for A = ±1.0 this anomalous peak 
disappear. 



B. Pair correlation function 

The nearest site correlation function between dimers 
can be obtained using a derivative of the free energy given 
by eg. pop . Then, in order to perform the derivative of 
free energy, we need to assume the parameter J' and J as 
independent parameter, similarly A' and A also are con- 
sidered as independent parameters. Despite for physical 
quantities, the term 2 J' + A' in the free energy be irrel- 
evant, this term cannot be neglected when we calculate 
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Figure 4: (Color online) Specific heat as a temperature de- 
pendency, for several values of anisotropic parameter. 



the correlation function as expressed below: 



«) = - 


1 dj 

2d J' 


x' dXo 
2Ao dx' 


[«)-- 


df 


, , z' 5Ao 
Ao dz' ' 


<ol)=- 


Idf 

2dJ 


X dXo 
2Ao dx ' 


yc<yl) = - 


df 
dA 


z dXo 
Ao dz 



(31) 
(32) 
(33) 
(34) 



Here, we present the relation of the correlation function 
in terms of the largest eigenvalue of the transfer matrix 
Ao. 



IV. THERMAL ENTANGLEMENT 

Another interesting property we consider in this 
work, will be the quantum entanglement of the Ising- 
Heisenberg orthogonal dimer plaquette model. As a mea- 
sure of entanglement for two arbitrary mixed states of 
dimers, we use the quantity called concurrence [31|, which 
is defined in terms of reduced density matrix p of two 
mixed states 



assuming 



C(p) = max{0, 2A,nax - tr^/R}, 



R = pcjy ® ayp*ay ® a^ 



(35) 



(36) 



where K^ax is the largest eigenvalue of the matrix y R 
and p* represent the complex conjugate of matrix p, with 
cr'' being the Pauli matrix. 

For the case of infinite chain, the reduced density op- 
erator elements[32,] could be expressed in terms of the 



Figure 5: (Color online) Concurrence for several values of 
anisotropic parameter as a function of T. Dashed line corre- 
sponds for a6-dimer while the solid line represents cd-dimer. 



correlation function between two entangled particles [33|, 
consequently the concurrence between a6-dimer becomes 



Cab = max{0, |(a>,-)| - -|1 + {a^oDW, (37) 
similarly the concurrence between cd-dimer reads 

C,a = max{0, |(a,VS)| - i|l + {a^aDW. (38) 

Surely, we can obtain also equivalent result using the 
approach described in reference [16J. 

It is worthy to mention that the zero temperature en- 
tanglement for a6-dimer is maximally entangled in TS 
and TT region Cab — 1, while, in all other region the 
a6-dimer becomes untangled {Cab = 0). Whereas the cd- 
dimer also is maximally entangled in same (TS and TT) 
region Ccd = 1; furthermore the DFA region now also be- 
comes an entangled region whose concurrence is given by 
Crd — I and it depends of J and Jo. It is worth to 

notice, when J = this region becomes untangled, which 
is perfectly coherent since the model reduces to orthog- 
onal plaquette Ising chain. Another special possibility 
is when Jo — 0, the orthogonal dimer plaquette reduces 
simply to independent entangled dimers. While the FM 
and AFM region are untangled Ccd = 0. 

In figure [5] is plotted the concurrence C as a func- 
tion of the temperature: for A = —1 the dashed and 
solid line corresponds to a6-dimer and cJ-dimer respec- 
tively. However, the a6-dimer entanglement vanishes 
for T K, 2.9 while, for cd-dimer the entanglement van- 
ishes at T sa 4.0. For weaker anisotropic coupling (A) 
this difference becomes more significantly, say i.e. for 
A — —0.24, the cd-dimer (solid line) entanglement van- 
ishes for T sa 2.25, whilst for a6-dimer the entanglement 
vanishes at T w 0.28 (dashed line). For a bit weaker 
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Figure 7: (Color online) Density plot of concurrence in units 
of Tth, darkest region corresponds to higiier entanglement, 
(a) Concurrence for afe-dimer as a dependency of J and A for 
a fixed value of Jo = 1. (b) Concurrence for cd-dimer as a 
dependency of J and A for a fixed value of Jo = 1. 



Figure 6: Threshold temperature Tth as a function of A, for 
fixed Jo = 1 and J = 1. Dashed line corresponds for a6-dimer, 
while the solid line represents cd-dimer. 



anisotropic parameter A = —0.22, the ab-dimer becomes 
untangled region for any temperature, while cd-dimer en- 
tanglement vanishes at T sa 2.24. 



A. Threshold temperature 



The threshold temperature can be obtain when ab- 
dimer concurrence becomes null (KcaCr^)! — 2^^ '^ 
(cr^cr^)l), and similarly for cd-dimer, the concurrence will 
become null when ((ctcCIJ)! = i|l -|- (cr^cr^)|). 

The threshold temperature is illustrated in figure [6] as 
a function of anisotropic parameter for fixed values of 
Jo = 1 and J = 1. The dashed line correspond for the ab- 
dimer threshold temperature, and in the low temperature 
limit, the threshold temperature leads to Tth = 2 — -v/S- 
Whereas the solid line represents the cd-dimer thresh- 
old temperature; in the low temperature limit it leads to 
Tth = a/5 — 2. It is worthy to highlight that the thresh- 
old temperature for cd-dimer exhibits a small reentrance 
around to A = 0.1. 

Another way to display the entanglement phase dia- 
gram is in units of threshold temperature as displayed 
in figure [71 the boundary between entangled region and 
untangled region is given by a solid (red) line. The black 
region correspond to maximally entangled region, and 
white region corresponds to untangled region, while light 
gray region corresponds to weak concurrence similarly 
dark gray corresponds to strong concurrence. The entan- 
gled phase diagram J against A in units of Tth, behaves 
quite similar for both a6-dimer and cd-dimer, although, 
for a6-dimer the untangled region appears already sig- 
nificantly for A ^ —3 and small values of J, while for 
cd-dimer the untangled region becomes relevant only for 
A ^ — 1 and small values of J. Certainly this is in agree- 
ment with the zero temperature phase diagram displayed 
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Figure 8: Density plot of Concurrence in units of Tth, darkest 
re^5[l|region corresponds to higher entanglement, (a) Concur- 
rence for afe-dimer as a dependency of Jo and A for a fixed 
value of J = 1. (b) Concurrence for cd-dimer as a dependency 
of Jo and A for a fixed value of J = 1. 



in figure [21a). 

Additional density plot of entangled phase diagram is 
displayed in figure [8| in terms of A and Jo , once again in 
units of Tth, the phase diagrams for a6-dimer is given in 
figEt^) I ^nd the phase diagram for cd-dimer is illustrated 
in fig[DJb) . Here, we observe the phase diagram between 
afe-dimer and cd-dimer are quite different, but still are 
closely related to the zero temperature phase diagram 
illustrated in figure [2l^b). Although the darkest (strongly 
entangled) region are very similar for both dimers. 



V. CONCLUSION 

In this work, we studied the dimer-plaquette Ising- 
Heisenberg chain, assembled between plaquette edge also 
known as orthogonal dimer plaquette pHlOl ■ Using the lo- 
cal gauge symmetry of this model, we are able to map 
onto a simple spin-1 like Ising and spin 1/2 Heisen- 
berg dimer model with single effective ion anisotropy. 
Thereafter, this model can be solved using the decora- 
tion transformation |25i427l . |29| and transfer matrix ap- 
proach. First we discuss the phase diagram at zero tem- 



perature of this model, where we found five ground states 
as illustrated in figure [5J one ferromagnetic, one anti- 
ferromagnetic, one tripled-triplet disordered and triplet- 
singlet disordered phase, beside a dimer ferromagnetic- 
antiferromagnetic phase. It is interesting to remark that, 
in the limit of pure Ising model this exhibits a frustrated 
region. Furthermore, are discussed the thermodynamic 
properties such as entropy, specific heat as well as the 
correlation function. For the specific heat curve we ob- 
serve an anomalous peaks close to the first order phase 
transition in the low temperature limit (figH]), for weak 
Heisenberg anisotropic coupling, but for strong Heisen- 
berg anisotropic coupling, this anomalous peak disap- 
pear. Additionally, using the nearest site correlation 
function it is possible to analyze the pairwise thermal 
entanglement for both afe-dimer and cd-dimer, as well as 
the threshold temperature of the entangled region are 
discussed as a function of the Hamiltonian parameters. 



There is some significant difference between both dimers 
(a6-dimers and cd-dimers) which is in agreement in the 
low temperature limit. However, for strong entangle- 
ment, both dimers are quite similar. As a consequence of 
this difference, we can illustrate one interesting result, re- 
garding to the reentrance type of threshold temperature 
for cd-dimer, despite for a6-dimer there is no reentrance 
temperature. 
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